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Many natural fibrous networks with fiber diameters much smaher than the average poresize can 
be described as three-dimensional (3D) random line networks. We consider here a 'Mikado' model 
for such systems, consisting of straight line segments of equal length, distributed homogeneously and 
isotropically in space. First, we derive analytically the probability density distribution p(r„o) for 
the 'nearest obstacle distance' r„o between a randomly chosen test point within the network pores 
and its closest neighboring point on a hne segment. Second, we show that in the limit where the 
line segments are much longer than the typical pore size, p(r„o) becomes a Rayleigh distribution. 
The single parameter a of this Rayleigh distribution represents the most probable nearest obstacle 
distance and can be expressed in terms of the total line length per unit volume. Finally, we show by 
numerical simulations that a differs only by a constant factor from the intuitive notion of average 
'pore size', defined by finding the maximum sphere that fits into each pore and then averaging over 
the radii of these spheres. 
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I. INTRODUCTION 

Many biological systems can be structurally described 
as random line networks. A typical example are gels 
that self-organize by the polymerization and subsequent 
crosslinking of filamentous proteins, such as collagen or 
fibrin. In order to characterize the stochastic geome- 
try of such systems, a frequently used parameter is the 
average poresize fpore- It is determined by finding, for 
a reprensentative fraction of network pores, the largest 
spheres that can be fit into that pores and then comput- 
ing the average of the radii of these maximum spheres. 
While this can be done numerically in a straight forward 
yet time consuming way, this definition of poresize is not 
suited very well for exact analytical calculations. There- 
fore, we suggest as an alternative measure the most prob- 
able nearest obstacle distance a for randomly chosen test 
points and show that it is directly related to the poresize. 

We start our investigations with a 'Mikado'-like net- 
work model that has two parameters, the length I of the 
hne segments and the volume density p of their centers. 
It is possible to compute the distribution p{rno) of nearest 
obstacle distances analytically in this model. In the limit 
of zero line length, the Mikado model contains the case 
of point networks. More interesting is the opposite limit, 
where / is much larger than the average poresize. In this 
case, the Mikado model converges towards a more gen- 
eral model that represents any random line network with 
a large persistence length. Indeed, the single system pa- 
rameter in this limiting case is the overall line density A, 
i.e. the total line length per unit volume. This parameter 
only sets the spatial scale of the network, and no other 
details matter for the distributions p(r„o) or W{rpore)- 
For example, a network composed of random circles with 
identical overall line density would yield the same uni- 
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FIG. 1. Collagen gel with a concentration of 1.2 mg/ml (Scale 
bar is 10 /im). Shown is the maximum intensity projection 
from a stack of 15 single confocal images recorded at a z- 
distance of 340 nm (total height = 5.1 //m). The fibers are 
straight on the scale of a typical pore size r^ore and their 
diameter is much smaller than rpore- Therefore, the structure 
of the system can be well approximated by a Mikado line 
network in the long fiber limit. 



versal Rayleigh distribution p(r„o) as the Mikado model, 
provided the radius of the circles is much larger than the 
average pore size. We compare these analytical results to 
a numerical simulation that is directly based on the ex- 
axt analytic geometry of points and lines, thus avoiding 
any possible artifacts arising from voxelation. 

After demonstrating perfect agreement of the simula- 
tions with the analytic results, we use the simulations to 
determine the poresize distribution for line networks of 
various density parameters A. As expected from scaling 
arguments, the average poresize Vpore = ca is simply pro- 
portional to (7, allowing us to determine the conversion 
factor as c wl.86. 



II. MODEL AND THEORY 

Distribution of nearest obstacles distances p(rno), 
accessible volume fraction Q{r) and pore sizes 

We consider random biphasic networks, in which ev- 
ery point of S-dimensional space either belongs to phase 
(pore, Uquid) or phase 1 (material, solid). In order to 
map out the stochastic geometry of the network, one can 
repeatedly choose a random point Rq = (x, y, z) within 
the 0-phase of the network and then find its 'nearest ob- 
stacle distance' rno{Ro), defined as the Euchdcan dis- 
tance from that point Rq to the closest point of the 1- 
phase (compare Fig. [2|a)). The network is then charac- 
terized by the distribution p(r„o) of the nearest obstacle 
distances. 

Closely related to p(rno) is the 'accessible volume frac- 
tion' Q{r), defined as the fraction of the 0-phase in which 
a sphere of radius r (from now on called a r-sphere) could 
be centered without overlapping the 1-phase (compare 
Fig. [2]jb)). In general, the dimensionless quantity Q{r) 
has the value Q{r — 0) — 1 and decreases monotonically 
for all radii r > 0. 

The complemental quantity 1 — Q{r) is the fraction 
of 0-phase for which an r-sphere overlaps the 1-phase. 
It corresponds to the probability that a random 0-phase 
point Rq has a nearest obstacle distance r^^ smaller than 
r, or 



Q[r) = Prob(rno < r) 
p{rno)drno- 



(1) 



The derivative of this equation with respect to r shows 
that Q{r) is just the negative cumulative probability of 
p(rno): 



p{r : 



dr 



Q{r). 



(2) 



While both quantities carry the same information about 
the network, the cumulative Q{r) is more convenient for 
analytical considerations, as will be demonstrated below. 

Another way to characterize pores of a network is to 
find the maximum sphere that fits to each pore and to 
define the 'pore size' Tpore as the radius of this maximum 



sphere. The concept is also illustrated in Fig. |2fc). 
denote the distribution of pore sizes by W{rpo 
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Random Line Networks: The Mikado model 

In the following we consider random networks in which 
the 1-phase consists of straight line segments of fixed 
length, with isotropic orientations and a homogeneous 
distribution throughout the 3D volume. We refer to this 
model as the Mikado model. 

Each individual line segment (LS) can be described by 
its center point and a unit direction vector. In order to 











^^^^ 


^^^^ 




FIG. 2. 2D illustration of various statistical measures 
used for networks of line segments, (a) Nearest obsta- 
cle distances rno (thin lines) for a few selected points (cir- 
cles), (b) Accessible volume (shaded areas) for spheres of a 
given radius, (c) Maximum spheres fitting into network pores, 
thereby defining the pore sizes r^ore- (d) A homogeneous, 
isotropic random distribution of straight line segments. The 
segments have a prescribed length I and their center points 
a spatial density of p. (e) Classification of line segments in 
the 2D Mikado model. 1-group (squares): Centers within r- 
sphere. 2-group (full circles): Centers outside r-sphere, yet 
with chance of overlap. 3-group (empty circles): Remote seg- 
ments without chance of overlap, (f) 2D sketch of an r-sphere 
(green), a concentric spherical shell of radius R (gray) and a 
specific point (red) within this shell. From all line segments 
centered at the red point, only those can intersect the r-sphere 
with orientations falling into a cone of apex angle ui. 



avoid ambiguities, we require that all unit vectors have a 
positive z-component and thus 'point upwards' (compare 
Fig. |2][d)). The two parameters of the Mikado model are 
the length I of the LSs and the volume density p = y 
of their center points, where TV is the number of line 
segments within a volume V . 

Consider first the extreme case / — 0, where all LSs 
degenerate into their center points, and place a r-sphere 
randomly into the system. Note that the configuration 
of LS-centers throughout the volume is a spatial Poisson 
process with 'event rate' that is identical to the volume 
density p. On average, the r-sphere will contain a number 
of 



(3) 



LS-centers. The probability Q{r) that not a single LS- 
center lies within the r-sphere is given by the Poisson 
probability for fc = events, which is 



Qir) = Poisson {fc = Q,nav = nav,i^Q{r)} 



(4) 

Therefore, in the case of the random point network the 
accessible volume fraction is given by 



Qir) 



- pr 



(5) 
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We now turn back to the general case / > 0. As before, 
we can write 



Qir 



Ar) 



(6) 



In order to compute Uavir), we note that with respect 
to a given r-sphere, the LSs can be classified into 3 groups 
(compare Fig. [2](e)): 

• 1-group with LS-centers inside the r-sphere. 

• 2-group with LS-centers outside the r-sphere, but 
yet with a possibility of intersecting the r-sphere. 

• 3-group with LS-centers too far away to touch the 
r-sphere. 

Only the groups 1 and 2 contribute to n^vir). The con- 
tribution of the 1-group is identical to the case of point 
networks above: 



(7) 



The 2-group consists of LSs with centers in a sphere 
of radius r + (1/2) around the center of the r-sphere. We 
now consider in more detail the ones in an infinitesimal 
spherical shell of radius R around the center of the r- 
sphere, with < i? < r + {l/2). This R-shell contains a 
number of 



dN' = pdV = pAnR^dR 



(8) 



candidates for intersection. Among them, only those LSs 
will actually overlap the r-sphere that have orientations 
within a certain cone (compare Fig. ^f)). This cone has 
an apex angle oiio ~ 2 arcsin(r/i?) and the corresponding 
solid angle is 

n{R) = 47rsin2(w/4) 



(9) 



47r sin arcsin (r/i?) 
2TT (l - Vl - {r/R)^) 



Since the total solid angle available for LS orienta- 
tions is fttot = 27r (according to our convention that 
all unit direction vectors are pointing upward), the in- 
tersecting LSs amount to a fraction of il{R)/iltot = 

(l - v/l - {r/R)^'j ■ We conclude that the average num- 
ber of actual intersections from LSs within the R-shell 



IS 



dN(R) = dN' (l - - (r/Ry) 



= Anp (l - - {r/R)A R^dR. (10) 



The total contribution from all LSs of the 2-group is ob- 
tained by integration over the relevant R-shells: 



l-R=r+{l/2) 

UavAr) = / dN{R). 
J R=r 



(11) 



This integral can be performed analytically. Using the 
abbreviation 



one obtains 



/(l+2;)-/(l) 



(12) 



(13) 



By adding the contributions of both relevant groups, 
"av(f) = nav,i{r) + Uavai'r), and using Q{r) = e-""-(''), 
we arrive at an analytic expression for the accessible vol- 
ume fraction in the Mikado model. Defining another use- 
ful abbreviation 



g{x) := 3 



/(I 



1 
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the result can be cast into the form 



(14) 



(15) 



It correctly contains the limit of point networks, since 
g{l/r) — ?> for I 0. All the differences between point 
and LS networks are included in the 'perturbation func- 
tion' g{l/r). 

From the accessible volume fraction Q{r), we immedi- 
ately obtain the distribution of nearest obstacle distances 
vifno — 'r) = ^^Qi''') in the Mikado model. With in- 
creasing Tno, this distribution starts with p{rno = 0) =0, 
develops a single peak and then decays exponentially for 
distances much larger than the average pore size Rpore 
of the network. 



Mikado model in the long fiber limit 

We next consider the case I ^ r, where the LSs are 
much longer than the typical distances of interest. Since 
pifno) is exponentially small for distances beyond the 
average pore size, this limit can also be interpreted as 
I 3> RpoTc- Note that this is a typical situation for net- 
works of semi-flexible fibers, such as collagen. 

It is straight-forward to show that in this limit the per- 
turbation function diverges as g(l/r) 
fore obtains 



1^. One there- 



Qi^rir) = e 



(16) 



which is the 'right half of a Gaussian bell curve with 
standard deviation 



(17) 



cr = l/^2TTpl 



The corresponding distribution of nearest obstacle dis- 
tances is a Rayleigh distribution 



Pi»r(r) 



-e 2 



(18) 
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FIG. 3. (a) Distribution of nearest obstacle distances in 3D 
networks of straight line segments, for three different density 
parameters A. Analytical predictions of the Mikado model in 
the long fiber limit (dashed lines) are compared to numerical 
simulations (solid lines). The unit of length was set equal 
to the linear size L of the simulation box, which in turn was 
equal to the length I of the line segments, (b) Distribution 
of nearest obstacle distances p{r) (line) and pore size distri- 
bution W{r) (line with symbols) in a 3D network of straight 
line segments, for a density parameter of X = WO/L^ . (c) Most 
probable obstacle distance a (squares) and average pore size 
'"pore (circles) as a function of the density parameter A. In the 
long fiber limit, the ratio is constant with rpore/<J ~ 1.86. 



The most probable nearest obstacle distance, i.e. the 
value of r at which p{r) is maximum, is given by a. We 
note that the accessible volume fraction in Eq.(16) de- 
pends only on the ratio r/a. 



Therefore, all nearest ob- 
stacle distance distributions pi:^r(j') should collapse onto 
a universal distribution when the distance r is measured 
in units of a. In the long fiber limit, a dense and a di- 
lute Mikado network cannot be distinguished from each 
other, if the spatial scale is unknown. 



Relating a to line density 

It is remarkable that in the long fiber limit of the 
Mikado model, the properties of the network are com- 
pletely determined by the parameter combination pi, 
which appears in the quantity a = 1/ y^2npl. 

Remembering the definition of p as the volume density 
of LS centers, we can write 



pl^yl=^=: A, 



(19) 



where Lfot is the total length of all LSs. The new density 
parameter A corresponds to the total 'fiber' length per 
unit volume. It follows that 



(20) 



Numerical test of the Mikado model 



In order to test the predictions of the Mikado model, 
we have simulated random line networks and compared 
the resulting numerical pir^o) with the analytical results 
above. 

In the simulation, each line segment (of constant length 
I) was numerically represented by its center coordinates 
and a unit direction vector, as depicted schematically in 
Fig. ^d). Initially, a list of N such line objects was 
generated, with the center points distributed randomly 
throughout a cubic simulation box of linear dimension 
L (with homogeneous density p — ^ = ^) and with 
random, isotropic direction vectors [? ]. 

The distribution p^r^o) was determined by randomly 
choosing K — 10^ test points Rk=i...K within the simu- 
lation box, finding the nearest obstacle distance r,io(^fc) 
for each test point and then computing a histogram of 
these distances. The distance r^o{Rk) is found by first 
computing the distances rkn between test point Rk and 
all the lines n of the network and then finding the small- 
est of those values. Note that the distance rkn between a 
point and a line segment can be obtained exactly (with- 
out any 'voxelization' required). 

For the numerical test of the Mikado model in the long 
fiber limit, we prescribed the density parameter A, set 
L = I = I and computed the required number of fibers 
as = We found an excellent agreement between 

the analytical prediction and the simulation (compare 
Fig. ^a)). 



Relation between the most probable nearest 
obstacle distance a and the average pore size Tpore 

For any concrete network, it is possible to compute the 
nearest obstacle distance rno(^o) for each spatial point 
Rq, resulting in a so-called 'Euclidean distance map' 
(EDM). The pore centers of the network can then be 
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defined as the positions Rq = Rmax of tlie local maxima 
of the EDM and the pore size distribution W{rpore) is 
the distribution of the distance values rp}re — rno{Rmax) 
taken at these local maxima. 

Based on our numerically exact simulation of random 
line networks, as described in Sectjllj we have computed 
the pore size statistics W{rpore) and compared it to the 
corresponding distribution p(rno) of nearest obstacle dis- 
tances. As expected, W{rpore) is peaked at a larger value 
than p(7'„o) (compare Fig. [Sj^c)). In the long fiber limit, 
the ratio Tpore/o' between the average pore size and the 
most probable obstacle distance is a constant, i.e. inde- 
pendent from the density parameter A of the network. 
This follows from the fact that the distribution p(rno) is 
universal in length units of cr. To demonstrate the con- 
stant ratio, we have plotted fporeW and cr(A) double- 
logarithmically (compare Fig. |3jd)). 



ically homogeneous distributions of straight line segments 
in 3D space. In the limiting case when the line segments 



are much longer than the average poresize 



the dis- 



tances of random test points to the nearest line seg- 
ment are distributed according to a Rayleigh distribution 



Pirno) 



(21) 



The most probable distance a (peak position of the dis- 
tribution) is determined by the overall line density A, i.e. 
the total line length per unit volume, by 



1 



(22) 



The average poresize Vpore, defined via the radii of max- 
imum spheres fitting into the pores, is proportional to cr. 



with ' 



pore 



1.86 cr. 



III. SUMMARY 
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